Search for: All records

We present an energy-specific Bethe–Salpeter equation (BSE) implementation for efficient core and valence optical spectrum calculations. In the energy-specific BSE, high-lying excitation energies are obtained by constructing trial vectors and expanding the subspace targeting excitation energies above the predefined energy threshold in the Davidson algorithm. To calculate optical spectra over a wide energy range, energy-specific BSE can be applied to multiple consecutive small energy windows, where trial vectors for each subsequent energy window are made orthogonal to the subspace of preceding windows to accelerate the convergence of the Davidson algorithm. For seven small molecules, energy-specific BSE combined with G0W0 provides small errors around 0.8 eV for absolute and relative K-edge excitation energies when starting from a hybrid PBEh solution with 45% exact exchange. We further showcase the computational efficiency of this approach by simulating the N 1s K-edge excitation spectrum of the porphine molecule and the valence optical spectrum of silicon nanoclusters involving 6000 excited states using G0W0-BSE. This work expands the applicability of the GW-BSE formalism for investigating high-energy excited states of large systems. 

Double excitations are crucial to understanding numerous chemical, physical, and biological processes, but accurately predicting them remains a challenge. In this work, we explore the particle–particle random phase approximation (ppRPA) as an efficient and accurate approach for computing double excitation energies. We benchmark ppRPA using various exchange-correlation functionals for 21 molecular systems and two point defect systems. Our results show that ppRPA with functionals containing appropriate amounts of exact exchange provides accuracy comparable to high-level wave function methods such as CCSDT and CASPT2, with significantly reduced computational cost. Furthermore, we demonstrate the use of ppRPA starting from an excited (N − 2)-electron state calculated by ΔSCF for the first time, as well as its application to double excitations in bulk periodic systems. These findings suggest that ppRPA is a promising tool for the efficient calculation of double and partial double excitation energies in both molecular and bulk systems. 

block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.